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The renormalization of quasiparticle (QP) dispersion in bilayer high- Tc cuprates is investigated 
theoretically by examining respectively the interactions of the QP with spin fluctuations (SF) and 
phonons. It is illustrated that both interactions are able to give rise to a kink in the dispersion around 
the antinodes (near (7r,0)). However, remarkable differences between the two cases are found for 
the peak/dip/hump structure in the lineshape, the QP weight, and the interlayer coupling effect on 
the kink, which are suggested to serve as a discriminance to single out the dominant interaction in 
the superconducting state. A comparison to recent photoemission experiments shows clearly that 
the coupling to the spin resonance is dominant for the QP around antinodes in bilayer systems. 

PACS numbers: 



The elucidation of many-body interactions in high- 
Tc superconductors (HTSC) is considered as an essen- 
tial step toward an insightful understanding of their 
superconductivity. Angle-resolved photoemission spec- 
troscopy (ARPES) has provided a powerful way to probe 
the coupling of charge QPs to other QPs or collective 
modes. Recent ARPES experiments unveil several in- 
triguing features in the dispersion, the QP weight, and 
the lineshape: (i) A kink in the dispersion was observed 
in both the nodal and antinodal regions 53, S S 0- 
(ii) The QP weight around the antinodal region de- 
creases rapidly with the reduction of dopings [a], while 
changes a little around the nodal direction [g. (iii) Af- 
ter disentangling the bilayer splitting effect, an intrinsic 
peak/dip/hump (PDH) structure was seen around the 
antinodal region both in the bonding (BB) and antibond- 
ing (AB) band |2,I3- (iv) The kink around the antinodal 
region seemly shows a different momentum, temperature 
and doping dependence from that around the nodal di- 
rection (v) The antinodal kink is likely absent (or 
very weak) in the single-layered Bi2201 These fea- 
tures, especially (i), imply that the QP is coupled to a 
collective mode. So far, two collective modes of 41 meV 
spin resonance 0,11 and ~ 36 meV phonon %^ have 
been suggested, but which one is a key factor res pon sible 
for the kink is still much debatable Q, [3, S S H M ^ . 
So, it raises an important question as if the electronic 
interaction alone is responsible for the intriguing QP dis- 
persion, or to what extent the antinodal (nodal) QP's 
properties are determined by the strong electronic inter- 
action. 

In this Letter, we not only answer the above impor- 
tant question but also present a coherent understanding 
on the above features by studying in detail the respective 
effects of the fermion-SF interaction and fermion-phonon 
interaction on the QP dispersion based on the slave- 
boson theory of the bilayer t — t'~J model. We find that 
though both couplings are able to give rise to the kink 
structure near the antinodal region in the QP dispersion. 



they differ remarkably in the following aspects, (a) The 
lineshape arising from the spin resonance coupling ex- 
hibits a clear PDH structure, while the phonon coupling 
would lead to a reversed PDH structure, namely the peak 
is in a larger binding energy than the hump, (b) The for- 
mer coupling causes a rapid drop of the QP weight near 
the antinodal region with underdoping, but the latter 
has only a very weak effect. Moreover, the correspond- 
ing coupling constant for the fermion-SF interaction is 
reasonable consistent with ARPES data, in contrast to 
the much smaller value for the fermion-phonon interac- 
tion, (c) The bilayer coupling plays a positive role in the 
occurrence of the kink in the case of the fermion-SF inter- 
action, but has a negative effect on the formation of the 
kink for the fermion-phonon interaction. These results 
suggest that the SF coupling be a dominant interaction 
involved in the antinode-to-antinode scattering. 

We will consider separately the interactions of fermions 
with the SF and phonons. Let us start with the bilayer 
t — t' — J model with the AF interaction included, 
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where a — 1,2 denotes the layer index, ta.a' = t, Ja.a' = 
J if a — a', otherwise, ta,a' = tp/2,Ja.a' = J±/2 
and i — j. Other symbols are standard. In the 
slave-boson representation, the electron operators Cja- 
are written as cja = b^jfja, where fermions fia carry 
spin and bosons bi represent the charge. Using the 
mean field parameters Xij = J2a < f}afj<y >== Xa, 
^ij =< Mfjl - filfj^ >= ±Ao, and setting b ^ VS 
with S the doping density (boson condense), we can de- 
couple the Hamiltonian (1). Its Fourier transformation is 

given by, = Ek.. e^fLf ft^ ^^ifif f^-fi + 
h-c.) + EkJ'5t±ke^'=^^/^f /^'j + h.c] + £„, with = 
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—2{6t + J'xo){coskx + cosky) — iSt' cos cos ky — /i/, 
Afc = 2J'Ao(cosfc^ - coskjleo = 4iVJ'(xo + A^). 
ti_k = tp{coskx — coskyY /A jlJI and J' — 3 J/8. Diago- 
nalizing the Hamiltonian, we get the AB and BB bands 
with the dispersion e'^^'^'^ = Ck ± 5t±^^. Then, the bare 
normal (abnormal) Green's functions of fcrmions Qi^'^^ 

{Gw^'^''), and the bare spin susceptibility x"'"'(aja' = 
A, B) are obtained. The physical spin susceptibility 
is given by x = Xo (^os'^i<lzc/2) + Xo sin^(<7zc/2), with 
Xo = X^^ + X"""" and Xq = X^"" + x""^- 

The slave-boson mean-field theory underestimates the 
AF correlation |0, Q , so we need to go beyond it and 
include the effect of SF through the random-phase ap- 
proximation (RPA), in which the renormalized spin sus- 
ceptibility is, 

X±(q, OJ) = Xo* (q, + (a^q ± ^±)Xo (Q, ^)/2] (2) 

However, in the ordinary RPA (a = 1), the AF corre- 
lation is overestimated as indicated by a larger critical 
doping density b ~ 0.22 for the AF instability than the 
experimental data (5c = 0.2 ~ 0.5 Thus, we use the 
renormalized RPA in which the parameter a is deter- 
mined by setting the AF instability at the experimental 
value be- The fermionic self-energy coming from the SF 
coupling is given by, 

S(^^^)(k) = ±_l^^[(Jq + J^)2^+(q)5(;t;'')(k-q) 
q 

+ (Jq- J±)V(q)t;f.'^Hk-q)], (3) 

where, the +(— ) sign is for the normal (abnormal) 
self-energy Ss(uj), and the symbol q represents an ab- 
breviation of q, iijJra- The renormalized Green's func- 
tion is G(^'^)(k,zw) = [(Gi-^'^^(k,icj))-i + (Ak + 
sl^'^))2Gi^'^)(k, with gI-^'^' = - e(^-^) - 

jn(-4,s)-|_i^ and the spectral function is obtained via 
A{k,J) = -(l/7r)Im[G(fc, w -I- i(5)]. To determine the 
QP weight z, we first write the Green's function as 
G(k, iw) — M/{iijj + Sq) + N/iibj ~ E^) and get z via 
z = M/[l + 5S]a/5w(wo)], with ujo the pole of G. The 
parameters we choose are t = 2J, t' — —0.7t, J = 120 
meV and J± = OA J Except in Fig. 4, the doping 

density is set at ^ = 0.2. 

Fig.l (a) and (b) display the calculated QP dispersion 
for the BB band obtained from the momentum distribu- 
tion curve (MDC) jl^ in the antinodal and nodal regions, 
respectively. One can see from Fig. 1(a) that the antin- 
odal kink appears in some range of parameters tp and a, 
such as tp = 2.0J and a = 0.43. The RPA correction 
factor a = 0.34,0.43 and 0.55 correspond to the criti- 
cal doping density of the AF instability Sc = 0.02,0.05 
and 0.09, while 6c — 0.02 ~ 0.05 is within the experi- 
mental range Meanwhile, the ARPES experiment 
indicated that t±^exp (corresponding to dtp, 6 0.2) is 
44± 5meV [i3, i.e., tp = 1.6 ~ 2.0 J here. Therefore, in 
the reasonable parameter range, the interaction between 
fermions and SF reproduces well the observed antinodal 
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FIG. 1: The MDC dispersion of fermions due to the inter- 
action with spin fluctuations (a),(b) and (c), and phonons 
(d). Figs. (a), (c) and (d) are the results at (k^jTv) (antinodal 
region), and Fig.(b) at {kx,kx) (nodal direction). The scat- 
tered points in Fig.(c) are derived from the EDC [T^ . 



kink. In contrast, Fig. 1(b) shows that no kink is present 
in the nodal region. This can be understood from the 
conservation of the momentum in the scattering process, 
namely the nodal-to-nodal scattering involves a smaller 
transferred momentum than Q = (tt, tt) where the AF 
spin fluctuation peaks. Thus, we will focus mainly on the 
antinodal region in the following discussion. In Fig. 1(c) 
we replot the MDC dispersion together with the EDC 
derived dispersion . Near and below the region where 
the kink appears in the MDC dispersion, the EDC dis- 
persion breaks into a two-part structure, a peak and a 
hump. Therefore, the appearance of the antinodal kink 
has an intimate relation to the PDH structure in the EDC 
plot, being in agreement with experiments [T^ . 

In fact, the kink may also be expected if fermions 
are predominantly coupled to other collective modes. 
A hotly discussed mode responsible for the antinodal 
kink is the out-of-plane out-of-phase O bucking Big 
phonon HII^- To take into account this kind of mode, 
we may have the total Hamiltonian by including the fol- 
lowing interaction in the slave-boson mean field Hamil- 
tonian Hm, 

^^P^^H q)/l[,./k+q,.(4 + rf-q) (4) 

k,q,a- 

where S and d are the creation and annihilation opera- 
tors for phonons, g(k, q) = go[ ^^{k)^x{k-q) cos(g^ /2)- 
$y(k)$j^(k - q) cos(q:c/2)]/A/cos2 g^/2 + cos^ qy/2 and 
the detailed form of can be found in Ref. |9l[l0|. 

We note that the vertex glk, q = 0) ~ cos(kx) — cos(fcy) 
and vanishes for all k at q = (tt, tt) |^, llO|. Thus, the 
fermions near (tt, 0) are strongly scattered by this inter- 
action. The corresponding fermionic self-energy is given 
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by, 

E,,„(k) = ±J-^|5(k-q,q)|2 
q 

i)o(q)[5i:t](k-q)+CKk-q)]. (5) 

where the Green's function of phonon is Do{q) = 
2L0ql\{iiJ)^ — (ijJg)^], and a dispersionless optical phonon 
(Big) wiU be taken as fiujq = 36meV j^, [i3- Unhke 
the fermion-SF interaction, the couphng constant here 
is not available now. We have tried various values and 
found that the well-established kink can be obtained 
when gQ « 0.2J ~ 0.4J if = 1.7J as shown in Fig. 1(d). 
This result is quite similar to that for the SF in Fig. 1(a). 
So, a mere reproduction of the QP kink is not sufficient 
to single out the main cause of the renormalization, and 
we must resort to other features. 

A meaningful difference between the effects of these 
two interactions on the dispersion can be seen clearly 
from a comparison of Fig. 1 (a) and (d) , namely, the inter- 
layer coupling has opposite impacts on the antinodal 
kink. It enhances the kink feature for the fermion-SF 
interaction; as shown in Fig. 1(a), when tp decreases to 
1.7J, no antinodal kink is observed even for the strongest 
AF coupling a = 0.55 which is in fact beyond the exper- 
imentally acceptable value. In contrast, the kink feature 
is weakened by the interlayer coupling for the fermion- 
phonon interaction; as seen in Fig. 1(d), the kink disap- 
pears when tp increases to 2 J from 1.7J with go = 0.2 J. 
Recent ARPES experiments revealed that a more pro- 
nounced kink is present in the multilayered BiSrCaCuO, 
in sharp contrast to the case in the single-layered one Q , 
which can be considered as an indication to favor the 
fermion-SF interaction in the bilayer system. Because the 
spin susceptibility in the bilayer system involves scatter- 
ings between layers, the self-energy [Eq.(3)] has a feature 
that the fermions in the BB (AB) band is scattered into 
AB (BB) band via the SF in the odd channel x~ ■ The 
so-called spin resonance (a sharp peak in Imx) appears 
around Q both in the odd and even channels 
(see Fig. 2(a)). However, it is more prominent in the odd 
channel, which is in agreement with very recent exper- 
iments [13, mainly due to the larger vertex aJq — Jj^ 
( Jq = —2 J) [Eq.(2)], compared to aJq -I- J_l in the even 
channel. On the other hand, the AB band (and its as- 
sociated flat band near Q) is much close to the Fermi 
surface compared to the BB band, as shown in Fig. 2(b). 
As a result, the fermionic self-energy for BB band is large, 
and consequently the renormalization is strong. We have 
indeed observed that the AB band is much less affacted 
by this coupling, so no kink is present in this band (not 
shown here). As increasing tp, the splitting between the 
AB and BB band pushes the flat portion of the AB band 
to be more close to the Fermi level [Fig. 2(b) and (c)], so 
enhances the scattering l2^. However, for the coupling 
to phonons, the AB and BB bands contribute to the self- 
energy in the same way[Eq.(5)]. In this case, though 
the increase of t„ increases the contribution from the AB 
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FIG. 2: (a) \mx{(\,uj) vs u at {ti,tt). (b) and (c) show the 
bare dispersion of the normal state quasiparticle for tp = 2J 
and 1.0 J, respectively. 
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FIG. 3; The lineshape of fermions at different k points. 
Fig. (a) is obtained when fermions is coupled to spin fluctua- 
tions. Fig.(b) shows the result arising from the coupling to 
the B\g phonons. The inset in Fig. (a) shows the lineshape 
for the bonding and antibonding bands at (0, tt) , separately. 
The inset in Fig.(b) is the lineshape for the bonding band at 
k = (0,7r) and (0.27r,7r). 



band, the decrease from the BB band overcompensates 
that increase, and thus the self-energy decreases with tp 
on the whole. 

In Fig. 3, we show the lineshape for different k points 
from (0, tt) to (0,0. 27r). For the fermion-SF interaction, 
both the BB and AB spectra near (0, tt) consist of a low 
energy peak, followed by a hump, and then a dip in be- 
tween, though the intensity of the AB hump is much 
weaker than its peak intensity [inset of Fig. 3(a)]; both 
develop their own PDH structure near (0, tt). When 
moving from (0, tt) to (0.27r, tt), one will see that the 
intensity of the BB peak increases, while the AB peak 
decreases rapidly. Eventually, only the BB peak can be 
seen near (0.27r, tt). These features are in good agreement 
with ARPES experiments 0|. However, the lineshape 
caused by the phonon-coupling displays a striking con- 
trast to those shown in Fig. 3(a) and in experiments ^3 
, namely, the peak is far below the dip and the dip is 
below the hump [Fig. 3(b)]. This is because the renor- 
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FIG. 4: The quasiparticle weight z of fermions and its cou- 
phng constant A at the antinodes arising from the couphng 
to spin fluctuations (a) and (b), and to phonons (c) and (d). 
The insets show those at the nodes. The scattered points are 
experimental data from Ref. Q and Ref. 0. 

malization to the QP peak in the BB band due to the 
phonons is rather weak, so the peak is almost unchanged 
comparing to the bare one. When moving from (0, tt) to 
(0.27r, tt), the QP moves to be near the Fermi level, an 
ordinary PDH structure is recovered [inset of Fig. 3(b)] 
because the hump arising from the phonon coupling does 
not change in the process. 

Figure 4 presents the QP weight z of fermions and the 
coupling constant A. Notice that the weight of the phys- 
ical electron is bz due to the condensation of holons in 
the slave-boson approach [l^ . For the fermion-SF inter- 
action, the weight decreases rapidly with underdoping, 
from nearly 0.42 and 0.39 at doping b = 0.24 (the bare 
value is 0.5) to be below 0.075 and 0.15 at (5 = 0.08 
for the BB and AB band, respectively. On the other 
hand, the weight along the nodal direction decreases very 



slowly, and it is still 0.7 even at (5 = 0.08. This ex- 
hibits the highly anisotropic interaction between fermions 
and SF and is well consistent with what is inferred from 
ARPES '3| and a recent argument based on the analysis 
of experimental data [23 . Moreover, the coupling con- 
stant obtained using z = 1/(1 + A) shows a reasonable 
fit to experimental data (3l| [Fig. 4(b)], while that at the 
nodal direction is much smaller than the experimental 
data [inset of Fig. 4(b)]. This consistence is signifi- 
cant, because we use the well established parameters in 
the t — t' — J model with only one adjustable parameter 
a being fixed by fitting to experiments. In contrast, the 
weight decreases much slowly for the fermion-phonon in- 
teraction as shown in Fig. 4(c), and the couphng constant 
A is nearly 3 ~ 5 times smaller than experimental values. 

Finally, we wish to make two additional remarks. First, 
because the spin resonance contributes little to the node- 
to-node scattering, a kink structure in the nodal direc- 
tion should be caused by the other mode coupling, such 
as air in- plane Cu-0 breathing phonon 1]. The dif- 
ferent momentum, temperature and doping dependence 
between the nodal and antinodal kink [feature (iv)] sup- 
ports this point of view. Second, since the interlayer 
coupling plays opposite roles in the antinodal kink re- 
spectively for the fermion-SF and fermion-phonon inter- 
actions, it is expected that, even though a rather weak 
fermion-phonon coupling may be present and lead to a 
weak antinodal kink-like behavior in the single-layered 
cuprates, as possibly seen in the experiment for Bi2201 
4], the coupling is too weak to affect significantly the 
lineshape which is mainly determined by the SF and was 
shown to exhibit the peak/dip/hump structure in the 
single-layered case [T^ ]. 

We are grateful to F.C. Zhang for many useful discus- 
sions. The work was supported by the NSFC (10474032, 
10021001, 10429401, 10334090), the RGC of Hong Kong 
(HKU 7050/03P), the URC fund of HKU, and partly by 
RFDP (20030284008). 



[1] A. Lanzara et ai, Nature 412, 510 (2001). 

[2] S. V. Borisenko et ai, Phys. Rev. Lett. 90, 207001 (2003); 

A. D. Gromko et ai, Phys. Rev. B 68, 174520 (2003). 
[3] T.K. Kim et al, Phys. Rev. Lett. 91, 167002 (2003). 
[4] T. Sato et al, Phys. Rev. Lett. 91, 157003 (2003). 
[5] A. Damascelli, Z. Hussain and Z. X. Shen, Rev. Mod. 

Phys. 75, 473 (2003). 
[6] P. D. Johnson et al, Phys. Rev. Lett. 87, 177007 (2001). 
[7] H. F. Fong et al, Phys. Rev. Lett. 75, 316 (1995). 
[8] P. Bourges et al. Science 288, 1234 (2000). 
[9] T. Cuk et al, Phys. Rev. Lett. 93, 117003 (2004). 
[10] T. P. Devereaux, T. Cuk, Z. X. Shen, and N. Nagaosa, 

Phys. Rev. Lett. 93, 117004 (2004). 
[11] M. Eschrig and M. R. Norman, Phys. Rev. Lett. 89, 

277005 (2002). 

[12] O. K. Anderson et al, J. Phys. Chem. Solids 56, 1573 
(1995). 



[13] P. A. Lee, Physica C 317-318, 194 (1999); J. Brinck- 
mann and P.A. Lee, Phys. Rev. Lett. 82, 2915 (1999). 

[14] J. X. Li, C. Y. Mou, and T. K. Lee, Phys. Rev. B 62, 
640 (2000). 

[15] A. P. Kampf, Phys. Rep. 249, 219 (1994). 

[16] Besides MDC, another intensity plot in ARPES is the 

energy distribution curve (EDC). Due to the abnormal 

lineshape of EDC MDC is more often used to analyze 

the dispersion. 
[17] D. L. Feng et al, Phys. Rev. Lett. 86, 5550 (2001). 
[18] M. R. Norman et al, Phys. Rev. B 64, 184508 (2001). 
[19] S. Pailhes et al, cond-mat/0308394 (2003). 
[20] Notice that the interlayer AF coupling Jx also enhances 

the coupling to spin resonance [see Eq.(3)]. 
[21] D. E. Sheehy, T. P. Davis, and M. Franz, 

cond-mat/0312529. (2003) . 



